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Abstract 

Recently, questions have been raised about the role of Lorentz invariance in false vacuum decay. 
It has been argued that infinities may arise in an integration over Lorentz-boosted final states. 
This suggestion motivates a Minkowski-space analysis of the decay rate. We attempt to illuminate 
'"pi features of the amplitude computation, and argue that the total rate including excitations is both 

finite and Lorentz invariant. 
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1 Introduction: The Prospect of Infinities in Vacuum Decay 



The problem of vacuum decay, both in the presence of gravity and without, has been studied for almost 
four decades. A general theory of vacuum decay in systems with large numbers of degrees of freedom was 
developed in [1 , 2] . Early efforts to calculate the false vacuum decay rate in field theory from a Minkowski 
space/WKB perspective appeared in [3], These calculations provided a rough picture, but they raised 
certain questions, including a question of finiteness; it was argued that it is necessary to integrate over 
Lorentz transformations of the basic process, leading to a factor of the volume of the Lorentz group in 
the total rate. The critical step was taken by Coleman (and Callan and Coleman) [4, 5], who formulated 
the problem at weak coupling in the language of the Euclidean path integral treated in the semiclassical 
approximation. The tunneling rate was given by the exponential of a large, semiclassical term (the 
WKB term) times a prefactor, obtained from fluctuations about a solution of the Euclidean equations 
of motion. The question of finiteness, it was argued, was solved; indeed, all infinities appearing in the 
calculation of the prefactor could be absorbed in standard field theory counterterms. In later work, 
Vilenkin and Vachashpati discussed the problem of constructing a wave functional describing vacuum 
decay [6]. They argued that their resulting construction was Lorentz invariant. 

Coleman and De Luccia subsequently developed the theory of vacuum decay in general relativity [7] . 
One of their remarkable observations concerned the decay of Minkowski space to a negative cosmological 
constant configuration; they found that such decays either do not occur, or lead to gravitational collapse. 
Most subsequent researchers have assumed that this indicates that the fate of Minkowski space, at least 
in the cases of apparent collapse, cannot be investigated within quantum field theory (see, e.g., [8]). An 
interesting attempt to model the process taking advantage of AdS/CFT duality appears in [9, 10]. 

Recently, Dvali [11] has revisited these questions and argued that decay of a false vacuum in field 
theory on Minkowski space does not make sense. Subsequently, Ref. [ ] delineated theories which might 
be discarded based on these considerations. 

The basic argument, like the old discussion of [3], has nothing to do with general relativity, but 
instead represents a new version of the argument that one is required to sum over Lorentz transformations 
of some decay configuration. Ref [11] notes that if a state can decay to a lower energy state through 
bubble formation, it can decay to a slightly larger bubble, with larger negative energy, and some positive 
energy particle. Such a configuration will not be Lorentz invariant, and it is argued that the total rate 
will be infinite after integrating over all Lorentz transformations of the configuration. 

One objection to this argument immediately comes to mind. Because the false vacuum is invariant 
under spatial translations but not time translations (it decays), it cannot be invariant under Lorentz 
boosts. Therefore the rate to nucleate a bubble along with a particle or excitation with classical total 
four-momenta (±2?, ±p) is not the same as the rate to produce a boosted bubble plus excitation with 
appropriately boosted four-momenta. This fact was mentioned in a footnote in [11], but not elaborated 1 . 

Ref. [11] speaks also of nucleating two bubbles, with momenta (±m, 0). This problem is of higher order in the 
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The issue can also be formulated in position space (or most precisely by considering the particle as a 
localized wave packet.) There we expect that the amplitude for production of the bubble and particle 
depends on multiple length scales, including the relative coordinate, the size of the bubble, and the 
lifetime: for example, if the particle is produced kiloparsecs away, one would expect the amplitude to 
be highly suppressed. Of course, these general arguments do not automatically preclude a divergence, 
only one proportional to the volume of the Lorentz group. Therefore we are motivated to study the 
production of excitations more carefully. 

As is well-known, the calculation of first quantum corrections to the imaginary part of the false 
vacuum energy was first done by Callan and Coleman in the Euclidean formalism [ ] . In the Euclidean 
computation, there is no signal of such an infinity at any order of perturbation theory about the bounce. 
Indeed, Coleman proved [ ] that the lowest action Euclidean configuration is always 0(4) symmetric. 
As a result, there is no need (at least for these bubbles) on the Euclidean side even for an integration 
over 0(4). On the Minkowski side, this translates into a classically 0(3, 1) (i.e. Lorentz invariant) 
nucleated state and its perturbations. If other non-0(4) symmetric bounce solutions exist and were to 
be included in the computation, because of the compactness of 0(4), the Euclidean computation would 
remain finite. For there to be an infinity, the Euclidean calculation would have to fail to capture the 
Minkowski space physics in a drastic way. 

Although the Euclidean computation is clean and straightforward, we can attempt to see what 
happens in Minkowski space directly by summing over the rates to nucleate all excitations about some 
basic bubble configuration. It would be very surprising if the Minkowski and Euclidean calculations 
differed in the final result. Certainly an infinite result would involve a loss of unitarity, or a failure of 
the continuation from Minkowski to Euclidean space. 

In an infinite Minkowski universe, in order to make sense of the false vacuum, we need to think 
of it as the limit of a sequence of long-lived states in which the field expectation value in increasingly 
large spatial regions lies at the local minimum of the potential. We can imagine preparing the universe 
in a "false vacuum" state from this sequence by heating and cooling the universe so that it transitions 
thermally to this state. We will assume that such 'boundary' issues do not have a strong impact on 
the lifetime of the approximate false vacuum state, and take as our starting point Coleman's original 
description of the Minkowski space tunneling probability. He noted that the analytic continuation of the 
0(4) symmetric bounce solution describes the nucleation of a perfectly spherical bubble at a particular 
space-time point (in Minkowski space). One can consider the nucleation of the bubble, along with an 
infinite possible set of excitations. We will focus principally on the thin wall approximation, which 
already incorporate the basic issues. We will describe a suitable set of basis states for the production of 
these excitations upon nucleation. The treatment of low energy excitations, as we will see, is complicated. 
The notion of the nucleation time, for example, is not extremely sharp 2 . But as our interest is in the 

semiclassical analysis, and requires an analysis distinctly more involved than we will consider in this paper, but, as we will 
see, there is no reason to expect an infinity in this case either. 

2 To quote Coleman, "Like all descriptions of quantum- mechanical processes in the language of the old quantum theory, 
this one must be taken with a large grain of salt; it will certainly lead us astray if we try to use it to describe measurements 
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role of highly excited states, this will not be an obstacle. 



Our goals, then, will be to: 

1. Define the basis states for our computation of the nucleation of excitations of the spherical bubble 

2. Develop a methodology for computing the nucleation of these states and determine if the rate is 



3. Understand the action of Lorentz transformations on the states and the total rate 

The rate to produce a symmetric bubble with an excitation has precisely the features of the parti- 
cle/bubble pairs of Ref. [11]. In particular, in order to produce a positive-energy excitation, it is necessary 
that the bubble, at nucleation, be larger than the critical bubble (corresponding to an excitation of the 
I = angular mode, which is of negative energy.) For the lowest lying excitations, the computation 
is challenging, but for the question of infinities, it is enough to consider high energy excitations, and 
for these the analysis is reasonably straightforward. We will see that, precisely because the system 
transits further in the effective potential for the bubble radius, these high energy excitations are cut off 
exponentially, with a decay rate T suppressed by 



where a is an 0(1) constant and AE is the energy of the excitation. At a classical level, the action of 
Lorentz transformations map bubble states with average radius greater than that of the critical bubble 
to new bubbles with even larger radius. Lorentz invariance of the sum is immediate, as can be seen 
by considering the problem from a "passive" viewpoint. Because the undisturbed bubble (no radial or 
positive energy excitations) is Lorentz invariant, it has a well-defined nucleation time in any frame, and 
takes the same form in any frame; as do its spectrum of fluctuations. As a result, the rate is the same in 
any frame. From the active viewpoint, the states simply furnish a representation of the Lorentz group, 
and no further states are required. The total rate is then expected to be finite, whether integrating over 
boosts or summing over excitations. 

Having understood the symmetric case, we will then discuss 0(4) non-invariant solutions of the 
Euclidean equations with false vacuum-false vacuum boundary conditions at t — > ±oo. To our knowledge, 
no one has exhibited such solutions. In fact, we will show that there is no 0(4) non-invariant solution 
in any situation where the thin wall approximation is valid (to the best of our knowledge this is a 
new result.) On the other hand, there is no reason to think that, more generally, such solutions don't 
exist. For any would-be 0(4) non-invariant solution (which also must have an odd number of negative- 
eigenvalue fluctuations in order to contribute to the decay rate), the puzzle is most stark. On the 
Euclidean side, one must perform an 0(4) collective coordinate integral, but this gives a finite factor of 
the volume of 0(4). On the other hand, the Minkowski space continuation of the solution (if it exists) 

made just outside the potential barrier. Nevertheless, it is very useful as an asymptotic description." [14] 
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no longer possesses classical 0(3, 1) invariance. An infinite set of distinct solutions can be constructed 
by acting with Lorentz transformations on a given solution, and we might think that we must somehow 
"integrate over all of these paths." But, as in the 0(4) symmetric problem, we should be precise about 
what to compute. If one considers the probability of materialization of the asymmetric bubble on 
some space-like surface, the situation is complicated. For a typical bubble, there is no time where the 
fields are static, for example, and determining a suitable complete set of states is not simple, let alone 
computing the amplitudes to nucleate each. Absent a clearly formulated computation, one cannot argue 
that there is a paradox. Indeed, as argued above in the limit where the bubbles are treated as particles, 
even making a correspondence between a family of classical Minkowski solutions and a complete set of 
quantum-mechanical final states, we would not expect to find equal tunneling rates. Furthermore, based 
on our experience with the symmetric bubble, which already takes into account non-Lorentz invariant 
configurations (just at one higher order in perturbation theory), we see no reason to think that the 
Euclidean computation should fail. 

The rest of this paper is organized as follows. In sections 2 and 3, we discuss the nature of the 
Minkowski problem as a WKB problem. (This subject has a long history, including [3, 1, 15, 6]). We 
focus on the symmetric bounce, recalling how the leading tunneling exponent is obtained, as well as 
the nature of the fluctuations. In section 4, we discuss the production of excitations in an adiabatic 
approximation. We first explain how such corrections arise in a conventional Born-Oppenheimer bound 
state calculation, and then in the problem of decay. We write explicit formulae for the amplitudes 
for decay with excitations. In section 5, we apply the insights from this analysis to the vacuum decay 
problem. We readily see that amplitudes for highly excited states are exponentially suppressed, and that 
the overall rate is finite. In section 5.3, we explain why these results demonstrate the Lorentz invariance 
of the final rate. In section 6, we consider 0(4) non-invariant bounces. Finally, in an appendix, we 
comment on a second fashion in which large Lorentz boosts may alter the decay process by correlating 
multiple bubble nucleations. 

Our work has overlap with that of [ ] . These authors also note that the consequences of Lorentz 
invariance are not obvious. But we differ in a number of details; e.g., they argue that "a perturbative 
calculation of the probability for nucleating a negative energy bubble accompanied by a single positive 
energy excitation ('blob') is divergent" but say that this is an indication of the breakdown of perturbation 
theory. The difference arises, in part, because of our focus on space-like surfaces; indeed, they argue that 
it is integration over different times that is the source of the problem. They also assert that there are 
no 0(4) non-invariant bubbles, and thus do not consider the puzzle raised in this case. They do provide 
an illuminating, explicit example of a number of the issues through the Schwinger process. 
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2 Setting of the Minkowski Problem 



The argument of Ref. [ ] applies perfectly well to the the thin wall limit, and the problem is most 
readily addressed there, as there is a natural separation of time and length scales. Following Coleman, 
we will consider a theory of a single real scalar field, with a quartic potential in zeroth approximation, 
and a small splitting (e) between the two degenerate states. In the Euclidean calculation, one finds an 
0(4) symmetric solution. The bounce solution changes from false to true vacuum at r ~ Rq, where the 
variation of the field occurs on a scale of order the inverse mass of the meson. In the thin wall 
limit, Rq is large, with size proportional to 1/e. The cross section of the wall is then that of a domain 
wall connecting the two, approximately degenerate, vacua. Si is the tension of this domain wall and is 
proportional to /i 3 , and in terms of it the critical radius is given by 

Ro= 3 -^. (2) 

The continuation of the Euclidean bubble 4>(r 2 + t 2 E — Rq) to Minkowski space is a Lorentz invariant 
evolution <\> c (r 2 — t 2 — Rq). The bubble is nucleated at t = with zero velocity (0 = 0). At the beginning, 
the bubble moves quite slowly; it takes a time of order Rq for it to double in size. So it is natural to 
expand the classical solution around t = 0: 

t 2 ■■ 

cj) c i(x, t) = <j> cl (x, 0) + — (f> d {x, 0) + . . . . (3) 

The equation for small fluctuations of <p at early times is given by 

d 2 

-^8(j>-V 2 5(i> + U{r)8<i) = Q , (4) 

where 

U(r)=V"(cl> cl (r,t = 0)). (5) 

5<t> can be expanded in eigenmodes of the differential operator 

H = -V 2 8<p + U (r)S(j) (6) 

Since the variation of U occurs principally at Rq, one can treat this problem analogously to familiar 
problems in molecular physics, separating variables and approximating the orbital angular momentum 
term by . Taking 

ty=-, (7) 
r 

we can start with the eigenfunctions of the radial operator 

d 2 



<}i , , V"{<l> d ) ) Xn = < X n. (8) 



This operator has an eigenfunction of negative eigenvalue 

d 



Xo(r,n) = -j^4> cl (r-R ,t = 0), uj 2 =-2/R 2 . (9) 
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Figure 1: Schematic form of the potential U(r). 



Xo can be combined with the I = 1 angular modes to make translations, which are zero modes of H . 
Higher modes of the radial operator will have eigenvalues of order the meson mass-squared and there is 
a continuum extending to oo (see Fig. 1 for the schematic form of U(r).) 

The modes of excitation can then be divided into two classes; those with frequencies comparable to 
Rq 1 and those with frequencies large compared to Rq . The slow modes (at t = 0) are those associated 
with the approximate translation zero mode of the solution: 

<p(S,t) = <p cl (r-R(t),0)+ J2 aoAm{t)Y em {Q)-xo(r - R) + ^t,m(t)Y tm {n) - Xn {r - R) (10) 

e>0,m n>0,e,m 

+ J dw 2 a u ^ )m (t)Y im (Q)^Xu(r - R) 

(The n = £ = m = mode is omitted from the sum; this is accounted for by the mode R(t). The 
continuum eigenvalues are labeled by lu 2 .) 

R(t) is a collective coordinate which describes the overall growth of the bubble. Near t — 0, the 
Lagrangian for this mode is 

A<jr 

L(R, R) = —S!R 2 R 2 - V(R) (11) 

where 

V(R) = A-kSxR 2 ~ t-ireR 3 (12) 
o 

For t ~ R, the motion is relativistic, and this Lagrangian is inadequate. We will return to the case of 
large (Minkowski or Euclidean) time, when the motion is relativistic, in section 3.1. 



7 



At small times and for large £ [l ^> 1), the a n g m obey the equation 



Q"nlm + UJ nlm a nlm — 



(13) 



where 



nim 




(14) 



For ^> (j^-) 2 , the period of the motion is short compared to the time required for substantial 

changes in R, and the system can be treated adiabatically for short times. These features of the classical 
problem suggest an adiabatic approach to high excitations in the quantum problem. 

Our interest is to consider tunneling with excitations of the various modes. We expect significant 
suppression with £, m and the size of the corresponding occupation numbers. To understand this problem 
we proceed in steps. 

1. We review how, for the dynamics of the mode described by R{t), a conventional WKB treatment 
reproduces the leading contribution to the tunneling rate found by Coleman in the thin wall 
approximation. 

2. We estimate the rate for production of excited states in the decay. For energies large compared 
to 1/Ro, we argue that this problem can be treated in a a Born-Oppenheimer approximation. We 
discuss an analog in a quantum problem with multiple degrees of freedom. 

3. From the experience with the analog problem, we obtain the dependence of the amplitudes on the 
energy of excitations in the field theory. We demonstrate that the amplitudes die off exponentially 
for energies of order [x or larger, and that the resulting tunneling amplitude is finite. 

4. We consider the question of Lorentz transformations and invariance. In a semiclassical approx- 
imation we can study the action of Lorentz transformations on bubbles larger than the critical 
radius (nucleated along with a positive energy particle) and see that the result is a new bubble 
configuration (along with a boosted particle.) The new bubble is larger, so the partial decay rate is 
suppressed compared to the rate into the unboosted state. Overall Lorentz invariance is maintained 
in the sum over all bubbles and excitations. 

In the next sections, we consider these issues in turn. 

3 Leading Tunneling Exponent 

For small t, it is easy to see that the solutions of the equations of motion obtained from Eq. (11) sweep 
out the solutions of the full Euclidean field theory problem. The bubble is nucleated where the potential 
vanishes, i.e. at Rq given in Eq. (2). Expanding the potential about this point, 



V*{R-Ro)V' 



(15) 
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where 



V = -SiRo. (16) 

So for small times (t <C R), 

R(t) = R + ^t 2 . (17) 
Let's compare with the behavior of Coleman's bubble at nucleation. This is 

MVi*-t*-Ro)*Mr-Jir)- ( 18 ) 

ZHq 

So for small times, we see that we describe the classical solution. More generally, we require a relativistic 
generalization. This is also the case for small fluctuations. 

3.1 Relativistic Treatment 

The n = modes correspond to displacements of the bubble wall: 

R(t, Q) = R(t) + 8R(t, n) = R(t) + a lm Y em (n). (19) 

tm 

5<f>= d 7^ L Yae m Ye m (n)- (20) 

tm 

An approximate treatment of wall evolution, in the classical case, has been developed in [17]. This 
approach allows a straightforward relativistic treatment. For the equation of motion for R, one finds 

R=^{l-R 2 f /2 -\{l-R 2 ). (21) 
For a zero energy bubble, this gives: 



{1 - R2)1/2 = 3 ^- (22) 

For small (Minkowski) times, this agrees with our earlier solution. The associated action is that of the 
spherical bubble. For large times, the motion becomes relativistic. For the Euclidean problem, \R\ grows 
at small R, behaving as 

1*1 = ^- (23) 

For the fluctuations, calling Ag m — ao£ m , one obtains the equation 

A(m + S (, ± (1 - #,V. - i) A„„ + "- A2 "^ + "- 2) A„„ . 0. (24) 
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For small times (R s=s 0), this agrees with our earlier equations for a. For large £, the system can be 
treated adiabatically. This approximation certainly breaks down for Minkowski times for which R ~ 1; 
the first derivative terms in the equation of motion dominate. But for Euclidean times, when R is large, 
the second and third terms in the equation arc comparable for small £, and the third term dominates 
for large £. So one expects an adiabatic treatment to be reliable for large £, and not too misleading for 
small £. These statements will carry over to the quantum treatment below. 



4 Higher Excitations: Adiabatic Approximation 



The discussion above is suggestive of an adiabatic treatment of the problem of production of excitations 
during tunneling.. To gain insight into such processes, it will be helpful to first consider some examples 
from quantum mechanical systems with small numbers of degrees of freedom. As a warmup, we write a 
formula for the first corrections to the Born-Oppenheimer approximation. We will phrase the problem in 
the language of molecular physics, but the generalization to other systems is immediate. In the molecular 
physics problem, one has nuclei of mass Mj, and coordinates Ri (which will abbreviate simply as R), as 
well as electrons with coordinates x a (which we will abbreviate simply as x) and mass m. One defines 
the zeroth order problem by solving for the electron wave functions for fixed nuclear positions, ip{x a ;R). 
The corresponding energy eigenvalues define an effective potential for the nuclei. One then solves for 
the wave function of the nuclei, ^(R); the full wave function is a product: 

$ Nn (R, x) = *„ >jV (i?)Vn(a;; R). (25) 

These solutions, of course, are not exact. They are eigenf unctions of the Schrodinger operator only if one 
neglects terms in the Schrodinger equation. Thinking semiclassically, one is neglecting the speed 
of the nuclei compared to that of the electrons. We can obtain the wave function to the next order of 
approximation by writing, say, for the ground state: 

<S> (R,x) = t>o{R)Mx;R)+ c Nn * N (R)il> n (x;R) . (26) 

Plugging in the Schrodinger equation, and keeping terms with one derivative with respect to R acting 
on tpo yields 

Now consider a system with a particle described by a coordinate R and a metastable minimum in the 
potential. Suppose that R couples to an oscillator, x, with frequency oj(R), a slowly varying function 
of R, with u(R — > 0) = bJ\ and u(R — > oo) = uj . We want to ask the amplitude to tunnel from the 
lowest state in the well to states with excitations of the oscillator. This situation is analogous to the 
field theory problem - by tunneling further along the R potential, one can create oscillator excitations 
of higher energy (arbitrarily high energy in the field theory problem). On the other hand, wc might 
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Figure 2: Left: A contour plot of a potential exhibiting a coupling between a harmonic oscillator x and 
a coordinate R with a metastable minimum (at the origin). For large R, the potential reaches a stable, 
flat value along x — 0, and the harmonic oscillator decouples from R. Right: the i?-potential for x = 0. 

expect that there is a price for these excitations; in the WKB approximation, one would expect a larger 
exponent (action) for highly excited states, suppressing the amplitude. 

To investigate this question, it is helpful to consider it as a scattering problem. We take the potential 
to be 



where Vb(0) = 0, and Vb(oo) = — Vq. We suppose that ^ is everywhere small, so that an adiabatic 
picture is appropriate, and that it goes rapidly to zero in the region where Vq(R) < 0; this insures that 
transitions between different oscillator states occur only under the barrier (this will not be the case in 
the field theory problem; we will discuss this issue in the next section.) A sketch of the potential and 
its profile for x = are given in Fig. 2. 

We consider an initial state with total energy E . We can write the incoming wave function in the 
form: 



71=1 

Here ipn(x; R) is the harmonic oscillator wave function appropriate to the n th excitation with frequency 
w(R), i.e. 



V{R,x) = V (R) + -K(R)x 2 



(28) 



OO 




(29) 



,2 



ip n (x;R) 



x Q 1/2 H n (x/x )e 



(30) 
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where 

x = (uj(R)m)- 1/2 . (31) 

To the right of the barrier, *&E ,n is the solution of the Schrodinger equation for R with potential 

V n = V (R) + (n+^)u(R) (32) 

with energy E , obtained in the WKB approximation. Provided one can ignore |g, behaves, in the 
forbidden region, as 

*E ,n( R ) ~ exp (- dr^2M(V n (r) - E )^j (33) 

where R n is the turning point of the potential V n (R). If the excitation number, n, is not too large, this 
can be estimated by studying the integration region near Rq. Calling R n — i?o + SR, 

nui(Ro) 

MM w (34) 

it is simple to estimate the additional contribution to the WKB exponent arising from the integration 
region between i?o and Rq + SR. For R smaller than Rq, we can estimate the additional correction to 
the WKB exponent by considering R < R — AR, AR ~ few 6R), and nu) -C Vb; then one can expand 
the square root in powers of uj/Vq and again estimate the corrections to the WKB integral. This yields 
the estimate: 

* Eo , n (R) = *o(i?)exp -n / dr- V K J== e o (35) 

\ Jro-ar 2 v /y '(i? ) x (r - i? ) J 

where c is a number of order unity, and this formula is valid away from Rq. We will estimate the size 
of these terms when we turn to the field theory. But it should be noted at this stage that the second 
factor dominates, and quite generally, we expect that under the barrier, R wavefunctions in states with 
oscillator excitations of energy AE will be suppressed by a factor 

£F~ e v 6 l 1 (36) 

In the regions in which |^ is not extremely small, the ^E ,n cannot be reliably computed in the WKB 
approximation. The factor 5", however, still describes the leading suppression of the amplitude for the 
particle to reach the barrier while still in the excited state, and tunnel through. 

We are interested in resonances associated with the lowest state in the well (i.e. we want to take Eq 
close to this energy); the coefficient of the pole is proportional to the width of the state (the distance off of 
the real axis is similarly proportional to the inverse lifetime). At leading order in the Born-Oppcnheimcr 
approximation, the lowest state has overlap only with the ground state of the oscillator, and the usual 
analysis of the decay width corresponds to taking n = for the initial state. However, we would expect 
that at higher orders in the approximation there will be overlap between the system's energy eigenstates 
and the leading-order product states with nonzero oscillator excitations. The decay rate into oscillator 
excitations will then receive perturbative contributions from overlap factors and exponentials of the form 
(36). 
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5 Application to the Field Theory Tunneling Problem 



In the previous sections, we have considered tunneling in a system coupled to an oscillator. We have 
seen that the final states can include not only the oscillator in its ground state but in excited states as 
well. Tunneling directly to these excited states is suppressed by the longer distance the system must 
tunnel through the forbidden region. Our model was grossly simplified by the requirement that ^ 
vanish beyond the barrier. In more realistic models, because of the exponential suppression of the direct 
tunneling, the most favorable way to produce excited states in the far future is through tunneling to the 
low lying states, followed by "jumps" to the higher levels. 

The field theory problem possesses certain features which are different than those of our quantum 
mechanics example. The potential, Vo, tends to — oo as R — > oo, and the frequencies of the low lying 
modes still vary in the allowed region. As a result, after the nucleation of the bubble, excitations are 
continually produced. Clearly to obtain the total rate for production of a given state at time t, it is not 
sensible to simply sum over all of these different processes happening at different times. Conceptually 
the simplest approach is to sum over all possible excitations at the moment of bubble nucleation. The 
notion of "moment of bubble nucleation" is not completely sharp, but is adequate for our considerations. 

From the quantum mechanics example, we have seen that for states for which the energy is not 
too large, there is an exponential suppression with energy (Eq. (36)). For very high energies, we would 
expect that the suppression is even stronger. In the field theory problem, 

M = \SiR 2 ; V(R) = 47rS'ii? 2 - ^Trei? 3 . (37) 

We will principally be concerned with the potential near R = Rq, so it is useful to note that 

M^SiRl; V'~SiRo; u ~ ^ £ ^ +1 \ (38) 

Rq 

Therefore, for large £, the suppression factor behaves as 

J~ e -[AE 3 /Si] 1/2 . (39) 

At any finite order of perturbation theory about the bubble, one does not expect exponential growth of 
the density of states, so we expect the overall rate to be finite. 



5.1 Low-Lying Excitations 



It is interesting to consider some of the types of excitations which might be nucleated along with the 
bubble. In Section 2, we identified a complete set of states for the field theory at the moment of 
nucleation. These were labeled by occupation numbers for levels u>,£,m, where, up to terms of order 
l/i?Q, the lu were eigenvalues of the domain wall problem in two dimensions. In addition to the zero 
eigenvalue, there may be some number of discrete eigenvalues, and a continuum above that. Consider 
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first the discrete states at n = 0. In comparison with the quantum mechanics problem, the frequencies 
are: 



y/t{i + 1) - 2 

W£m = R • ( 4 °) 

The I = mode corresponds to the collective coordinate R; the £ = 1 modes correspond to translational 
collective coordinates, about which we will say more shortly. 

Overall, one has a suppression of the total rate controlled by the sum: 

f(£)e (S i H B) 1/2 . (41) 

where / varies more slowly with i than the exponential. This sum is convergent. The low-£ large-n 
terms in the sum cannot be treated in the Born-Oppenheimer approximation, but we do not expect a 
divergence from them. 

5.2 Particle Emission with the Bubble 

Now consider the continuum. Here we encounter directly the question of production of very high mo- 
mentum/energy particles. We again expect suppression, but we need to pay attention to questions of 
energy and momentum conservation. 

We first consider the problem from the point of view of the mode expansion. Emission of a particle 
of energy-momentum = (E,p) corresponds to emission of a state in the continuum; the energy- 
momentum must be compensated by either the emission of another continuum state, or of one of the 
bound excitations. Focusing on the bound modes, in linearized approximation, and at t = 0, 



d 3 xd^ cl doS ( p + O(<50 2 ). (42) 



This receives contributions only from the I = 1 mode. In this way, energy and momentum can be 
conserved in the emission of a high momentum particle by 

1. turning on a non-zero ao,i, m — 

2. varying the initial bubble radius, 

R = Ra + AR (43) 

with 

dV 

E =8R AR - (44) 

Note, in particular, that the kinetic energy of the bubble (the J d 3 x^(do(f>) 2 term in the energy) is 
suppressed (it is of order i^gz)- Momentum and energy can be conserved in other ways, for example, 
turning on other modes and working to higher order in the fluctuation parameters, a n j. m - 
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For the excitation of modes with 

n.r\ i n = 



vt 

ao,i,o = -5- ; Ai?^0, (45) 



we have 



l\ = j d 3 X d z <f> cl d 5^> = / '' ' c~ir - R) (46) 



AlTV * Z?2 f , ,,2 47r «» R 2 

— R / dr^o = -g-Si-Ro , 



£ = J^A-R = -4 7 rS'ii? Ai? . (47) 

It is straightforward to work out the Lorentz transformation of the fields through first order in p. For 
transformations along the z direction, there is a shift in the £ — 1 component and AR: 

AR 1 

P; AR^ AR - -j3-vR Q , (48) 
Ro 3 

corresponding to 

E -> E + /3P; P -> P + f3E, (49) 

as they must for a solution of the equations of motion. Note that for f3 antiparallel to v, Eq. 48 simplifies 
to v — > u(l + Ai? — > Ai?(l + pi). In addition, the boost generates an £ = 2 mode, but this contributes 
to the energy-momentum only at higher order in p. At higher orders, one must be careful to work out 
the solutions to the equations of motion to higher orders in t at small time. 

Another way to produce a high momentum particle and conserve 4-momentum is to consider the 
nucleation of the spherical bubble and particle at rest, and apply a large Lorentz boost to the system. 
If the stationary bubble is larger than the critical bubble, it is not Lorentz invariant, and its rest 
energy/invariant mass should be the negative of the particle mass. In Ref [11], it was pointed out that 
if the nucleation rate for the boosted system depends on the bubble-particle relative 4-momentum 
only as p 2 , then an integration over boosts will diverge. However, we expect that boosts will in general 
excite various modes of the bubble, which we can calculate in the classical limit using the classical time 
evolution of Eq. 21. In particular we expect boosts to lower the energy of the bubble further, which 
implies a growth of the bubble radius (£ = mode) as we boost. Then, from our previous analysis, 
we expect exponential suppression of the rates. Indeed, as stated in the introduction, the breaking of 
Lorentz invariance in the rates is expected since the false vacuum initial state is not Lorentz invariant. 

When the particle is at rest, the bubble must satisfy (to leading order in m/(SiR 2 l )) 

777 

Now define Ri = Rq + AqR and consider Lorentz boosts of the system. Using the small-time equation 
of motion, valid to order t 2 , we find that boosted bubble at t = satisfies to order P 2 

r 2 - (i~3p 2 -y 2 - ^ = ■ (51) 
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A R/tf o = 1 

i i i i 1 i 



0.0 0.2 0.4 0.6 0.8 



Figure 3: Growth of a bubble radius under Lorentz boosts of parameter j3. The unboosted bubble is 
spherical and static with radius Rq + AqR = 2Rq. 

Here we are using a membrane picture, which is equivalent to the field theory calculation above where 
we considered bound modes. The £ = mode of the configuration has indeed grown, with 



and the total AR = A R+ApR. We see that at order /3 2 there will be a contribution to the wavefunction 
suppression. To go to large boosts requires solving the full relativistic equation of motion for the spherical 
membrane, which can be done numerically. In Fig. 3 we plot the growth of ApR as a function of f3 for 
the case of an initial AR = i?o, and observe that the growth of R becomes rapid at large j3. Note that in 
the numerical analysis AR does not need to be small. Eventually the linear approximation used in the 
rate suppression will break down, but since the potential falls in only a polynomial way, the exponential 
behavior of the suppression will continue. 

5.2.1 Suppression of the Tunneling Amplitude for Large Momentum 

Given that AR > for emission of a particle with non-zero momentum, there is a contribution to the 
suppression of the tunneling amplitude of order 



A R 



^/3 2 A R , 



(52) 




(53) 
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where = \p\ if only £ = 1 is activated (as in the first example above), and /(|p|) is a more 

complicated function related to the function in Fig. 3 in the case where p is obtained by boosting both 
a particle at rest and a static spherical bubble. 

Because the initial false- vacuum state is not Lorentz invariant, this suppression is consistent with the 
relativistic invariance of the underlying theory. We have stressed that we are interested in the amplitudes 
for emission of bubbles with excitations at time zero. We have seen that bubbles are mapped by Lorentz 
transformations into bubbles of different radius, with some linear motion and some distortion. As we 
will explain in Section 5.3, one may also consider the effects of Lorentz transformations on the particle 
wave function. This wave function may be considered as a localized wave packet (amplitudes depend 
on the distance from the wall) , and the Lorentz transformation of the wave packet changes not only the 
average momentum but also the location and shape of the packet. 

5.2.2 Two Analog Models 

The field theory description of these solutions, found in the linearized approximation, and their energy, 
momentum, and boosts is straightforward, but it is interesting to provide descriptions in terms of small 
number of degrees of freedom. Two analogies suggest themselves. 

First, for AR > and no excitations of non-zero n or £, the bubble has the structure of a particle at 
rest with negative energy. Thus it is tempting to describe it a negative mass particle, with Lagrangian 

i 771 

L = -m A /l-P 2 ; E = — - ; m < 0. (54) 

V 1 — v 2 

Small boosts of the particle at rest yield a negative rest energy and a negative kinetic energy, ^m/3 2 . 
The corresponding energy and momentum can be balanced against the energy and momentum of any 
emitted particle. In this description, it is not immediately obvious whether this negative "kinetic energy" 
is associated with a suppression of tunneling. 

Alternatively, imitating the conventional collective coordinate procedure for solitons, we can intro- 
duce collective coordinates for translations, 

x — >• x — xq (t) (55) 

with corresponding lagrangian (ignoring, for the moment, mixing with the collective coordinate R(t): 

L = -AttS^x 2 - V(R) (56) 

i.e. a non-relativistic particle of mass 47rS'ii?Q. This interpretation is consistent with the expression for 
the momentum above; it is, of course, valid only for small velocities. 

In this latter description, Lorentz transformations induce both a change in = v, and in AR. 

Av = -3-^ (57) 

Hq 
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In particular, in the limit of small j3 (Galilean transformations), the transformation of the velocity is not 
that which one has for ordinary particles; the induced velocity is not the parameter of the transformation, 
and there is a change in the radius. This analogy has other limitations as well. In particular, the "mass" 
in the non-relativistic limit for the collective mode is not the invariant mass of the configuration. 

Both of these pictures are useful in understanding aspects of this problem, but each has its limitations 
and they should be used with some care; it is safest to work directly with the field variables. 

5.3 Lorentz Invariance 

Let us now comment on Lorentz invariance from two points of view. From the passive viewpoint, Lorentz 
invariance of the total rate is straightforward for the 0(3, 1) symmetric bounce. We have described a 
calculation of the tunneling rate in Minkowski space. The principle limitation of the calculation is the 
difficulty of the study of low I excitations, but we do not expect this to generate a new divergence. 
The calculation involved the study of fluctuations about a Lorentz invariant configuration. In any other 
frame, if we define t — as the time at which the time derivatives of the background configuration 
vanish, the calculation is the same. In other words, rather trivially, as anticipated by Coleman, the 
calculation is Lorentz invariant. 

From the active viewpoint, as, for example, used in Sec. 5.2, we need to Lorentz transform states 
at t = 0, our emission point. If we sum over a complete set of states, the Lorentz transformations will 
just reorganize the sum. But there is no reason that the amplitude to produce the transformed field 
configuration needs to be the same as the amplitude to produce the original configuration. Indeed, the 
boost requires knowledge of the time evolution of the states, which will certainly be complicated in the 
full quantum theory, but in the semiclassical limit we have seen that boosts generate £ = 0, 2 deformations 
at 0(P 2 ) when acting on larger-than-critical spherical bubbles. We have argued that because of the £ — 
deformations, individual rates are suppressed under boosts, even though the sum over rates is Lorentz 
invariant. 

We can also note some features of the action of boosts on particle wave packets that might be 
nucleated along with the bubbles. Consider a high-momentum wave packet (semiclassical) located at 
t = at a point on or close to the bubble wall. High-momentum particles are nearly free, so we can 
consider the effects of Lorentz transformations on the wave packet alone. For simplicity, we take the 
wave to move along the x axis and limit our considerations to Lorentz transformations along this axis. 
A wave packet with average momentum p$, group velocity v g = (jgjf) > an d centered on the wall at the 
time of bubble nucleation takes the form 




(58) 



The Lorentz transform of this configuration behaves at t — as 




(59) 
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For large boosts, this corresponds to a wave packet at the nucleation time centered at 

x = i? V / l - v 2 /{l - v g v), (60) 

i.e. well-inside the bubble for sufficiently large boosts. Therefore, the particle is nucleated further 
and further away from the increasingly large-radius bubble wall, and so we expect further amplitude 
suppression with increasing boost parameter (in addition to the previously-discussed suppression induced 
directly by the growing radial mode.) 

From our analysis of this section and the rest of Section 5, we conclude that there is no reason 
to expect that the Euclidean analysis is in some way misleading. While the Minkowski calculation is 
distinctly more challenging than its Euclidean counterpart, it also leads to results which are finite and 
Lorentz invariant. 

6 0(4) Non-Invariant Bounces 

It is interesting to consider the behavior of 0(4) non-invariant bounces. Coleman et al. proved that the 
lowest action solution is always 0(4) invariant. We will demonstrate below (Section 6.1) that all thin 
wall solutions are 0(4) invariant. To the best of our knowledge, non-0(4) invariant solutions have not 
been exhibited explicitly, but we know of no general argument that non-symmetric solutions do not exist. 
Such solutions continue to Minkowski space solutions which are not 0(3, 1) invariant, and so pose the 
puzzle of Ref. [ ] quite sharply: does one have to integrate separately over all 0(3, 1) transformations 
of the solution (and the small fluctuations about it), and does this yield an infinite result? On the 
Euclidean side, the integral over the collective coordinates associated with the broken 0(4) is manifestly 
finite. But in continuing to the Minkowski side (if such a continuation can be defined sensibly- some 
rotations of a non-0(4)-symmetric configuration will break t — > — t symmetry, leading to imaginary field 
values under naive analytic continuation) , one of the Euclidean angles becomes a boost parameter, and 
there would appear to be an infinity of configurations to which one can tunnel. 

In the case of 0(4) symmetric bubbles and (non-0(4) perturbations around them) we had to deal 
with the question of what precisely one should calculate in a Minkowski space analysis. We suggested 
choosing a particular space-like surface (the approximate moment of bubble nucleation) and summing 
over all possible configurations on that surface. We saw that the Lorentz transforms of such configurations 
were complicated. What one might describe as nucleation in one frame with emission of a high momentum 
particle on the wall would have a quite different description in another. This is certainly the case for a 
classical, non-0(4) invariant solution. Indeed, there would be no preferred moment of nucleation (in an 
arbitrary frame), even in the approximate sense we considered in the previous section. If one chooses an 
arbitrary space-like surface as the nucleation surface, just as in the case of high momentum excitations 
we have considered earlier, the amplitude for the appearance of a bubble and the Lorentz transformed 
bubble need not be the same on that surface, just as in the case of non-Lorentz invariant excitations of 
the symmetric bubble. Of course, lacking explicit solutions, it is hard to sharpen these questions. But, 
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as in the symmetric case, it would be quite shocking if the Euclidean computation were not to produce 
the correct amplitude. 

6.1 Only 0(4) Invariant Solutions in the Thin Wall Approximation 

In the thin wall approximation, we can show that the 0(4) invariant solution is the only solution that 
satisfies the equation of motion with the boundary condition that the field tends to the false vacuum at 
infinity. 

In this limit, the size of the bubble is much larger than the thickness of the wall. Therefore, very close 
to the wall, any solution will look like a domain wall solution. The bubble solution is then characterized 
by the shape of the wall. Let B be the interior of the bubble and dB its boundary. We consider only 
the case where dB is topologically a three sphere S 3 . The Euclidean action can be written as 



where Si is the tension of the wall, e the (negative) energy density of the false vacuum, and g the 
determinant of the induced metric of the flat 4 dimensional metric on dB. 

Note that the action is just the difference between the area of the bubble and the volume it encloses. 
It is well known [18] that a solution that extremizes this action has constant positive mean curvature 
and the sphere is the only embedded surface that satisfies this condition. Therefore, we conclude that 
the sphere S 3 is the only type of solution. The radius of the sphere can be determined by extremizing 
the action for the bubble of radius r. This singles out the Coleman's bubble as the unique solution. 

7 Conclusions 

There is a history of concern about Lorentz invariance in tunneling computations in Minkowski theories, 
and of divergences arising from Lorentz transformations of particular tunneling configurations. To 
address these issues, we explained why it is crucial to formulate the question of what one actually 
calculates. There is no analog of an S'-matrix for such problems, and the initial "state" , given that it 
lives only a finite time, is not itself Lorentz invariant. For the thin-walled bounce, we argued that it is 
natural to think in terms of a bubble nucleation time, and of the complete set of excitations which may 
be produced at this time. Quite generally, we demonstrated that production of excitations of energy AE 
directly through tunneling (as opposed to through excitation from low lying states after the tunneling 
event) is suppressed by at least 



(where S\ is the bubble wall tension). As a result, there is no divergence with high energy states, e.g. 
with the emission of very high momentum particles or production of high £ modes of the bubble. 




(61) 



e -(AB 3 /Si) 1/2 



(62) 
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We began with the assertion the equivalence of the Minkowski and Euclidean computations usually 
follows from simple unitarity and analyticity considerations. While our study has elucidated the structure 
of the Minkowski computation and the absence of infinities, the actual real-time calculation of the process 
is complicated, and we have only dealt with some aspects here. Indeed, one lesson, if it was not obvious 
already, is that the Euclidean computation is far simpler than the Minkowski one. 

Ultimately, it is necessary to consider general relativity in understanding questions of tunneling. 
Coleman and De Luccia showed there are situations where tunneling does not occur in gravitational 
theories. In particular, a zero or small positive c.c. state will not necessarily tunnel to a state with 
negative cosmological constant (crunch space-time). On the other hand, in a suitable limit, gravitational 
corrections are unimportant in determining the decay rate. In fiat space theories, we have seen no 
difficulties with the conventional analysis. Indeed, we have seen that the Minkowski version of the 
computation, while complicated in detail, has a structure which can be readily understood. 

Finally, we noted that 0(4) non- invariant bounces would be interesting to study from this perspec- 
tive. Our analysis of the excitations of the symmetric bubble, however, indicates that no fundamentally 
new difficulties are to be expected. Furthermore, we showed that such asymmetric bubbles do not exist 
in theories in which the thin wall approximation is applicable. 
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A Large Boosts and Past Bubbles 



As explained in the Introduction, we do not expect the partial decay rates into final states related by a 
Lorentz boost to be the same. The boost-non-invariance can be understood from the dependence of the 
rates on multiple length/time scales in the problem. In this appendix we consider one effect of the time 
scale associated with the intrinsic finite lifetime of the false vacuum. 

In the picture of Ref [11], we consider matrix elements of the form (3 0i is a boost generator) 

M = (dole^^V^*-* ^) , (63) 

where \d Q ) is a classically-non-0(3, l)-invariant final state and \i) is the false vacuum. We can rewrite 
M as 

M = (do\e- mt - to) e luJo * 3 °a \i) , (64) 

which tells us that we can either study the projection of the time evolution of the initial state onto 
various final states boosted around the event (t,0), or equivalently we can study the projection onto 
one final state of the time evolution of various boosts of the initial state around the event (to, 0) on the 
initial data surface. We will now analyze the second case. 

In particular, since we are interested in decay rates per unit volume, we can ask how large the boost 
has to be to move a bubble into the past lightcone of a given point in the future of the initial data 
surface. Once a bubble is present in the past lightcone of the point, the amplitude to nucleate a new 
bubble at that point goes to zero. Let us choose our initial surface (where the false vacuum is prepared) 
at t = and require that all points below the initial surface are free from bubbles. Our future point is 
the event P = (r, 0) ; this may be regarded as the farthest point in the future for which it is still probable 
not to have a bubble in the past lightcone between P and the initial data surface. Under a Lorentz 
boost about the origin, some events inside the past lightcone are moved out and are replaced by events 
that were previously outside the lightcone. A two-dimensional illustration is given in Fig. 4. 

If the 4-volume of points above t = that are pushed into the past lightcone by a Lorentz trans- 
formation is of the order of the 4-volume of the past lightcone of P above t = 0, then the boost may 
be considered the smallest transformation that will with non-negligible probability introduce a bubble 
somewhere into the past lightcone of P (including prior to t = 0.) After some geometry, we arrive at 
the equation 



7tt fv((5 - 2v)v + 11) 



, — 3 arcsinh (try) I (65) 

6 V (v -l) 2 (v + l) y "J y ' 

where the left-hand side is the 4-volume of the past lightcone and the right-hand side is the 4-volume 
of points that pass into the lightcone under boosts determined by velocities of magnitude less than or 
equal to v. The right-hand side is obtained by integrating the area of the shaded region in Fig. 4 over 
boosts in the other spatial directions. Since r is the only dimensionful parameter in this problem (in the 
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Figure 4: A spacetime diagram showing events (shaded region) that enter the past lightcone (solid green) 
of the future point P located at (1,0) under a series of Lorentz boosts about the point (0,0). Under 
one such boost the dashed blue lines are mapped to the solid green lightcone. The shaded region is 
bounded by the initial data surface (thin black line at t — 0) in order to show only those points which 
may contain a nucleation and that can be mapped somewhere into the past of P (including prior to the 
initial surface) under a boost of magnitude less than or equal to the boost mapping A to P. 

limit of a small, coincident bubble+particle final state), as long as it is finite and nonzero, it drops out 
of the equation. We can then solve for a constant value of v, independent of the lifetime of the system. 
We find a numerical value of v max ~ 0.17. 

Consequently, even if all quantum mechanical objections raised in the main text above could be 
overcome, the contribution to the lifetime from states obtained by boosts greater than v max is likely to 
be negligible. 
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